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Resumen 


El trabajo propone un modelo matemático para la administración óptima de la 
explotación de acuíferos llamado MADA. Dicho modelo está basado en la inclusión 
de un modelo de simulación del flujo de agua subterránea (AQUIMPE) dentro de 
un algoritmo de programación lineal. Se plantea la concepción teórica de MADA 
realizando una especial distinción a la inclusión de AQUIMPE dentro del modelo 
propuesto. Finalmente se emplea MADA en la evaluación de la explotación de la 
“Cuenca Sur de La Habana”, haciendo uso de un Plan de Uso del Agua (PUA), 
donde queda demostrado que la distribución de la explotación empleada para 
garantizar el PUA no evalúa las posibilidades de intrusión salina en el acuífero y a 
tal efecto se brinda una solución. Además, se evalúa la posibilidad del acuífero de 
satisfacer la demanda bajo un patrón anual de lluvia del 75% de probabilidad, en 
un periodo de cinco años secos consecutivos. 
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Introducción 


Desde la década de los años sesenta, los sistemas 
de agua subterránea han sido estudiados 
mediante modelos matemáticos de simulación, 
los cuales se han venido desarrollando de una 
forma vertiginosa. En los últimos años, los 
modelos de simulación han sido combinados 
con técnicas de optimización para determinar 
la mejor estrategia a seguir en la administración 
del recurso hídrico, con un objetivo específico 
y un conjunto de restricciones (Galloway et al., 
2003). 

Estas técnicas de optimización pueden 
ser clasificadas en dos grupos: las técnicas no 
determinísticas, dentro de la que se pueden 
mencionar a las redes neuronales y los algoritmos 
evolutivos fundamentalmente; y las técnicas 
determinísticas, dentro de las que se encuentran, 
entre otras, la programación lineal, la cuadrática 
y la programación dinámica (Marrero, 1985). 
Para resolver estos problemas de optimización, 
muchas veces se aplica de forma aislada alguno 


de estos algoritmos, sin embargo, cada vez es 
más frecuente que se realicen combinaciones 
entre ellos para evitar obtener una solución no 
deseada o para disminuir el tiempo de cómputo. 
Por ejemplo, Hsiao y Chang (2002) proponen un 
método novedoso para encontrar la distribu- 
ción de pozos de bombeo, y el caudal a extraer 
en cada pozo y en cada intervalo de tiempo, que 
garantice un mínimo de costo global del sistema, 
donde se incluyen los costos de operación o 
bombeo, y los costos fijos o de emplazamiento. 
Este método combina un algoritmo genético 
para la determinación de la cantidad y ubicación 
de los pozos de bombeo, y otro algoritmo 
llamado Programación Dinámica Diferencial 
con Restricciones (CDDP) para la determinación 
de los costos dinámicos o de operación. 

Otros ejemplos de aplicaciones realizadas en 
el mundo en la administración de los recursos 
hidráulicos pueden ser consultados en Gutiérrez 
(2008). 

En Cuba se produjo un avance muy importan- 
te en esta rama del conocimiento con la creación 
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del modelo matemático de simulación del flujo 
subterráneo AQUIMPE (Martínez, 1989), en el 
Centro de Investigaciones Hidráulicas (CIH). 
Numerosas han sido las aplicaciones de este 
modelo en acuíferos nacionales e incluso en el 
extranjero (Gutiérrez, 2008). 

Sin embargo, sólo con una herramienta de 
simulación no se pueden resolver con rigor 
científico un conjunto de problemas referentes 
a la administración de la explotación. 

En Cuba, una de las primeras aplicaciones 
de patrones óptimos de explotación prolonga- 
da en un acuífero fue planteado por Dilla 
(1993). En este trabajo se emplea el método de 
la matriz respuesta o funciones tecnológicas 
algebraicas (FTA), donde finalmente se utiliza 
un método de optimización, la Programación 
Lineal, y de esta forma se resuelven varios 
problemas de optimización de naturaleza 
lineal. 

Posteriormente, Mesa (2004) da solución 
al problema inverso de estimación de los 
parámetros hidrogeológicos, conductividad 
hidráulica (K) y coeficiente de almacenamiento 
(S), aplicando un algoritmo genético de 
convergencia global llamado Shuffled Complex 
Evolution (SCE) y fundamentado en el 
modelo AQUIMPE. En este trabajo también 
se recomienda emplear el método SCE en 
la resolución de problemas referentes a la 
administración de acuíferos. 

Gómez (2009) incluye el modelo SCE en 
la tecnología AQUIMPE, para la estimación 
o el ajuste automático de los parámetros 
hidrogeológicos del acuífero, K y S. 

En la presente contribución se discute el 
modelo de administración MADA, basado en la 
inclusión del modelo AQUIMPE dentro de un 
algoritmo de Programación Lineal, y se aplica 
al acuífero de la Cuenca Sur de La Habana. 


Fundamentos y estructura general del 
modelo AQUIMPE 


El modelo matemático AQUIMPE en esencia 
resuelve numéricamente la ecuación diferencial 
parcial de segundo orden, de tipo parabólica, 


que surge al analizar el flujo bidimensional, 
impermanente, lineal, y libre o confinado en 
medio poroso, con término fuente-sumidero y 
condiciones de contorno filtrante o impermea- 
ble (Martínez, 1989): 


d pr d pr oh ¿50% , 
Ox dx) dy dy ot 


Donde: 


carga hidráulica [L]. 

espesor saturado del acuífero [L]. 
conductividad hidráulica [LT"]. 

coeficiente de almacenamiento, adimen- 


DANS 


sional. 


= 


fuente de entrada o salida vertical al 
elemento, es decir, explotación o recarga. 


La solución de esta ecuación puede expre- 
sarse en la obtención de la carga hidráulica (o 
piezométrica), como variable dependiente de 
las coordenadas del espacio bidimensional (x, 
y) y del tiempo (+). 

AQUIMPE (1) 
mediante el método de los elementos finitos 


resuelve la ecuación 
(MEP), utilizando el triángulo cuadrático como 
elemento, y aplicando la aproximación de 
Galerkin (Martínez et al., 2000). 

Con AQUIMPE se pueden simular acuíferos 
confinados y libres a escala regional, lo que 
incluye el tratamiento del acuífero visto en 
planta, aunque también se pueden resolver 
problemas vistos en una sección vertical. 
Además permite hacerlo a escala local, para 
simular, por ejemplo, pozos de bombeo y 
campos de pozos. Asimismo, es posible simu- 
lar su interacción con un cuerpo de agua 
superficial, siempre que exista comunicación 
entre ellos. 

La estrategia utilizada para la modelación 
del acuífero es subdividir el área total de 
la formación acuífera en distintas zonas, 
atendiendo a las propiedades hidrogeológicas; 
posteriormente estas zonas son subdivididas 
en elementos triangulares, donde en cada uno 
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de éstos las propiedades hidrogeológicas son 
iguales. Al emplearse el triángulo cuadrático 
para la discretización espacial, la solución en 
cada elemento se obtiene en seis nodos y por 
lo tanto ésta puede ajustarse a una superficie 
cuadrática, aspecto que permite representar 
con menor número de elementos en la malla 
una región, logrando la misma precisión de 
los resultados que si se emplearan triángulos 
lineales (Martínez, 2006). Además, con el 
triángulo cuadrático se puede modelar mejor 
el medio cársico como elemento volumétrico 
representativo de un medio poroso equiva- 
lente. En este sentido, la estrategia de modela- 
ción es definir la hipótesis de sustitución 
del medio cársico real por un medio poroso 
equivalente, cuyas propiedades de comporta- 
miento se determinan definitivamente durante 
la calibración del modelo (Cabrera, 2007). 

No es objetivo plantear todo el desarrollo 
del MEF aplicado a la ecuación (1) (Martínez, 
1989). Por eso, se ilustra el sistema matricial 
expuesto en la ecuación (2), obtenido después 
de la resolución de las integrales que se gene- 
ran con la aproximación de Galerkin a la 
ecuación (1): 


ne 
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Donde: 
Ao é ; 
[R] = m, HK.[Kab |, ¿ se conoce como matriz 


de rigidez del elemento. 

A.: área del elemento e [L?]. 

K;: conductividad hidráulica del elemento e 
[LT"]. 

[K.,]: ¿: matriz simétrica cuadrada de orden 6, 
cuyos valores son los mismos para todos los 
tiempos de un elemento e, pero diferentes para 


todos los elementos. 


Se Ao 


180At 
At: cantidad de unidades de tiempo en el 


[El = [ Kcte]¿, matriz constante. 


intervalo de tiempo. 

matriz constante para todos los 
elementos en todos los tiempos. 

(H.Y..,: valor de las h de los nodos del triángulo 
o elemento e; cada término del vector está 
referido a uno de los seis nodos del triángulo y 
depende del tiempo k [L]. 

WM yJ%..: igual al anterior, pero en el tiempo 
previo [L]. 

Qs': caudal superficial, en forma de intensidad 
que entra al elemento e, recarga debida a la 
lluvia en el tiempo k [LT]. 

Qp; caudal puntual que entra o sale al nodo ¡ 
de explotación en el tiempo k [L*T”]. 

QI': caudal lateral que entra o sale del acuífero 
por el lado ¡ de un elemento e [L?T"]. 

L: longitud del lado ¡ [L]. 

(000111), ,: vector columna, los tres primeros 
valores son los referentes a los nodos principa- 
les y el resto a los nodos secundarios. 

[0...100) 


elementos como nodos de la malla y en él tiene 


| Vector columna que tiene tantos 
valor uno los referentes a los nodos donde hay 
gasto puntual. 
(114), ,: vector columna de tres elementos; 
garantiza que en los lados donde hay gasto 
lateral se reparte 1/6 de este caudal a cada 
nodo principal y los restantes 2/3 al nodo 
secundario. 
n: número total de nodos de la malla. 
nl: número de lados de triángulos con salidas o 
aportes laterales (zonas del acuífero que tienen 
frontera con un río, etcétera). 
np: número de nodos con gastos puntuales. 
ns: número de triángulos con entradas o salidas 
superficiales. 
ne: número total de elementos (triángulos). 
Jj: contador de nodos con gastos puntuales. 
1: contador de lados con gasto lateral. 
e: contador de elementos. 
: tiene un valor de 2/3 si se obtiene por el 
MEF. 

Cuando se aplica esta expresión a un 
acuífero confinado, el problema se restringe a 
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resolver k veces el sistema matricial dado en 
la ecuación (2), el cual se resume en el sistema 
de ecuaciones lineales mostrado en la ecuación 


(3): 


MO mt Fa = Ur (3) 


Donde [M] es una matriz cuadrada de 
orden n, en la que cada fila está asociada con 
uno de los nodos de la malla. Lo mismo sucede 
con las columnas. En el vector h se guardan 
las variables a calcular en cada tiempo y f es 
el vector de los términos independientes donde 
se almacena el efecto de los caudales extraídos 
e incorporados al modelo en cada nodo. 


Estructura del modelo de programación 
lineal 


El asistente matemático MATLAB (García et 
al., 2005) posee una función llamada Linprog, 
que resuelve el problema de programación 
lineal mediante el método Simplex. Esta 
función es usada en el modelo propuesto y 
en su forma más general presenta la siguiente 
estructura: 


[x, feval |] = linprog (c, A, r, Asgr logo lb, ub, x,) (4) 


Donde c, A y r son las matrices que in- 
tervienen directamente en el problema, x es el 
vector de variables de decisión que se obtiene 
y c es un vector de costo o de ponderación 
de cada una de las variables de decisión. 
Con la ecuación (4) se resuelve el siguiente 
problema: 


mín cx (5) 
Sujeto a restricciones de desigualdad: 
Ax<r (6) 
Restricciones de igualdad: 
Aa (7) 


Con fronteras inferiores y superiores para 
cada variable de decisión: 

Ib<x<ub (8) 

El vector x, representa una solución inicial 


del problema y en el término feval se almacena 
el valor evaluado de la función objetivo. 


Funciones objetivo y restricciones del 
modelo de administración 


En el modelo MADA se pueden tomar en 
cuenta un conjunto de funciones objetivo 
(Dilla, 1993). 

Antes de analizar las funciones objetivo 
se define el término npa como la cantidad de 
nodos de la malla numérica con gasto puntual 
a optimizar. 


a) Maximizar la suma de las cargas en los 
nodos con gasto puntual a optimizar npa 
del modelo en todos los tiempos: 


T npa 
F.O=máxY Y h' 


kl ja1 


(9) 


b) Minimizar la suma de los abatimientos si 
calculados por la diferencia entre los niveles 
en el estado inicial y al final del horizonte 
de planeamiento o tiempo final T en todos 
los npa: 


F.O= mín Y si (10) 
Pl 


De esta forma se satisface la demanda con la 
menor afectación al acuífero. 

Por otro lado, la tercera función objetivo 
concibe la máxima extracción posible. 


c) La mayor capacidad de extracción de 
un acuífero se obtiene al maximizar los 
caudales en los nodos de extracción a 
optimizar npa para todos los tiempos: 
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T npa 
F.O=máx Y Fo! (11) 
k=L Jal 
En general, las restricciones pueden 


resumirse en tres tipos: 


1. La suma de los caudales extraídos en cada 
periodo por todos los nodos de bombeo 
del modelo a optimizar npa debe igualar 
o superar una demanda prefijada en el 
periodo. Este conjunto de restricciones 
genera T ecuaciones como la (12): 


(12) 


2. El caudal de extracción de cada nodo de 
bombeo del modelo a optimizar npa no 
debe sobrepasar la capacidad máxima 
proyectada para los pozos asociados con el 
nodo; aunque al menos debe entregar una 
demanda local especificada. 


Dicha serie de restricciones genera T x npa 
ecuaciones (por cada límite), es decir, 2T x npa 
ecuaciones en total: 

Q¡mín< Oj <Q;máx (13) 

3. Los niveles piezométricos en cada nodo de 

bombeo del modelo a optimizar npa para 

cada tiempo deben ser controlados, esto 

sugiere que estén por encima de un valor 

límite prefijado por el modelador. En este 
caso se generan npa ecuaciones: 


hi> he (14) 


mínj 
Donde: 


hs nivel piezométrico en el nodo de extracción 
a optimizar j. 
h* 


mín 


E nivel crítico en el nodo de extracción a 
optimizar ¡. Normalmente es un dato, valor 
constante en el tiempo. 


Modelo de administración (MADA) 


Las funciones objetivo definidas en las ecua- 
ciones (9), (10) y (11) son lineales y en este 
sentido es posible emplear un algoritmo 
de programación lineal para resolver los 
problemas de optimización que sugieren di- 
chas expresiones. 

Elmodelo de administración MADA se basa 
entonces en introducir todas las ecuaciones 
del modelo de simulación AQUIMPE dentro 
de un algoritmo de programación lineal, 
empleando para tal efecto la función Linprog 
de Matlab. 

Es importante plantear que este modelo será 
efectivo para acuíferos confinados o libres que 
en su modelación puedan ser tratados como 
confinados (Gutiérrez, 2008). 

El modelo de programación lineal quedaría 
de la siguiente forma: 


e Función objetivo: cada vez que se ejecute 
el modelo se desea resolver un problema 
específico que tendrá que ver con algunas 
de las expresiones (9), (10) u (11). Además, 
se debe tener en cuenta que las variables 
de decisión del modelo son los caudales 
de explotación en los nodos que se quieran 
optimizar en cada intervalo de tiempo, Q' y 
también los niveles piezométricos en todos 
los nodos en cada intervalo de tiempo. 

En esta función objetivo se debe declarar 
un vector de costo c (Gutiérrez, 2008), el 
mismo tiene una fila y ((n x T) + (T x npa)) 
columnas. Todos sus elementos son cero, 
excepto los que participen en la función 
objetivo que se esté analizando, que toman 
un valor de —1. 


e. Restricciones de igualdad: éstas vienen 
expresadas en Linprog con la ecuación (7) 
y tienen una gran importancia, ya que es 
en este sistema de ecuaciones donde se 
incluyen todas las ecuaciones del modelo 
de simulación, es decir, el modelo no 
encuentra los niveles en los nodos tiempo a 
tiempo con la ecuación (3), como se realiza 
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con AQUIMPE en el problema directo, sino 
que tiempo a tiempo el sistema (3) se va 
incluyendo dentro de otro sistema mayor, 
que es el representado en la ecuación (7), 
quedando una gran matriz A, con (n x 
T) filas y ((n x T) + (T x npa)) columnas, 
donde están los coeficientes que afectan a 
todas las variables de decisión. El vector a 
queda con ((n x T) + (T x npa)) filas y una 
columna. 

Cuando se resuelve el sistema planteado 
con la programación lineal se obtienen 
todos los niveles piezométricos en los 
nodos y en todos los tiempos, de una vez. 


e. Restricciones de desigualdad: las restric- 
ciones de desigualdad son inecuaciones 
que permitirán introducir en el modelo las 
importantes restricciones que se abordaron 
en (12) (13) y (14). La restricción (12) 
se introduce mediante la matriz A y el 
vector r dados en la inecuación (6), la cual 
está concebida por la función Linprog 
de Matlab como una inecuación del tipo 
menor o igual que. Cuando las inecuaciones 
restrictivas de (12) son del tipo mayor o igual 
que, se multiplican todos los elementos de 
A y rpor-1. 


La matriz A queda con (T) filas y ((n x T) + 
(T x npa)) columnas y los vectores x y r quedan 
con ((n x T) + (T x npa)) filas y una columna, y 
(T) filas por una columna, respectivamente. 

Para introducir las otras restricciones, 
dadas en las ecuaciones (13) y (14), se utilizan 
los vectores de fronteras inferiores y superiores 
[Ib) y [ubj, cada uno de estos vectores tiene 
dimensiones de ((n x T) + (T x npa)) filas y 
una columna. El valor de cada elemento 
de los vectores es el mínimo o máximo, 
respectivamente, que puede tener su variable 
de decisión asociada. 


Discusión 
Las experiencias que se tenían en Cuba en la 


administración Óptima de la explotación del 
agua subterránea empleaban las funciones 


tecnológicas algebraicas (las cuales se obtienen 
mediante un trabajo previo del usuario 
con el modelo de simulación AQUIMPE), 
para posteriormente aplicar el algoritmo de 
programación lineal. Sin embargo, en MADA 
esto no es necesario, ya que en el proceso 
de optimización se tiene al modelo de flujo 
propiamente dicho, lo cual es una ventaja, ya 
que mejora todo el trabajo de modelación a 
realizar por el usuario, así como el tiempo de 
cómputo. El precio a pagar es la complejidad 
de la concepción e implementación de este 
modelo, que se le suele llamar modelo acopla- 
do (Pérez, 1995; Yaragal y Nagaraj, 2007). 

Otras de las virtudes del modelo MADA es 
que por emplear AQUIMPE, que usa el trián- 
gulo cuadrático como elemento finito, tiene una 
muy buena precisión en los resultados que se 
obtienen. Por lo general, los modelos basados 
en el método de los elementos finitos que se 
encuentran en la literatura no usan el triángulo 
cuadrático como elemento finito. 

Además permite: 


e Declarar recargas al acuífero en todos los 
tiempos y cada elemento. 

e Declarar flujos laterales de entrada o 
salida. 

e Declarar nodos con gastos puntuales 
conocidos (que no serán optimizados). 

e Declarar nodos con cargas fijas en todos los 
tiempos. 

e Declarar nodos con cargas variables en 
todos los tiempos. 

e Limitar inferior y superiormente cualquie- 
ra de las variables de decisión (cargas y 
caudales), e incluso a combinaciones de 
ellas. 

e Amplia variedad de 
incluidas animaciones. 


salidas gráficas, 


Modelo regional del acuífero “Cuenca 
Sur” de La Habana 


Características de la zona de ubicación del 
acuífero 


La zona de estudio se encuentra en el sur de 
la provincia de La Habana, abarcando un área 
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de 830 km? (ver figura 1). Existe en esta zona 
una gran carsificación de las rocas del complejo 
carbonatado de mioceno, una diversidad de 
las formas del carso, compleja relación con 
la tectónica, una abundante alimentación al 
acuífero por filtración y poco espesor de la 
cobertura de las rocas permeables de la zona 
de aireación, entre otros factores, que crean 
condiciones favorables para la formación de 
considerables recursos de agua subterránea 
(Llanusa, 1997). 


Las condiciones hidrogeológicas de la 
región son bastante complejas (Llanusa, 1997) 
debido a: 


1. Gran heterogeneidad y anisotropía de las 
rocas con diferentes grados de carsifica- 
ción. 

2. Cuenca abierta con descarga al mar. 

3. Compleja composición química del sistema 
acuífero desarrollado en las calizas y dolo- 
mías del mioceno. 


Y Número de triángulos 
e 


Y Nodos principales 
o 


Y Nodos secundarios 
o 


Y Dique Sur de La Habana 


nde > 


rd 


Figura 1. Ubicación de la zona de estudio y discretización del acuífero con la numeración de nodos y triángulos. 
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Guira de Melen 
A 


.- A 


Y Pozos de Y Nodos | Nodos 
O observación 2 principales 2 secundarios 


Y Dique sur Y Triángulos Y Niveles piezométricos 


N /V del agua subterránea 


Figura 2. Mapa digital de isolíneas de niveles del agua subterránea correspondiente al estado inicial. 


4. Extracción intensa e irregular del agua. Se 
realiza principalmente en el periodo seco. 


A lo largo de la costa y con una faja 
de ancho de aproximadamente uno o tres 
kilómetros, predomina el suelo cenagoso. El 
cultivo principal en la región es la caña, se 
observan también cultivos menores y pastos. 
Hacia el sur se tiene la vegetación típica costera 
representada por distintas especies: mangle, 
patabán y otras. 

Orográficamente la región es muy llana, las 
elevaciones oscilan entre 100 m al norte y 3-5 m 
al sur. La red fluvial no está bien definida, y es 
prácticamente imposible delimitar las cuencas 
hidrológicas. Hacia el sur existe una amplia 
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red de canales, cuya alimentación proviene 
en lo fundamental de las aguas subterráneas 
(Llanusa, 1997). 


Modelo conceptual propuesto 


Toda la modelación matemática del acuífero 
está basada en un conjunto de abstracciones en 
síntesis de las características físicas reales del 
sistema natural que será estudiado con fines 
de aplicación práctica. Los criterios del modelo 
conceptual se exponen a continuación: 


e El acuífero puede ser simulado con un 
modelo bidimensional no estacionario en 
planta, esto quiere decir que sus propieda- 
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des hidrogeológicas no varían en dirección 
vertical. 

e El medio filtrante puede simularse como si 
fuera un medio continuo poroso isotrópico 
y heterogéneo. 

e Enel modelo de administración, el acuífero 
será simulado como confinado. 


Discretización de la zona de estudio 


En la conformación de los elementos de la 
malla discretizada se tuvo en cuenta que tanto 
la geología como las propiedades hidrogeo- 
lógicas se mantuvieran homogéneas en esos 
elementos. En la zona paralela a la costa se 
emplearon elementos que cubrieran el ancho 
de esta franja cenagosa, no pudiéndose utilizar 
elementos más pequeños por la carencia 
de información. La zona de estudio quedó 
dividida en 69 triángulos, con un total de 158 
nodos, como se observa en la figura 1. 


Condiciones de frontera 
La zona objeto de estudio se encuentra en la 


vertiente sur de la provincia de La Habana. Por 
el norte limita con el parteaguas del sur de la 


Cuenca de la Laguna de Ariguanabo y parte 
de la divisoria central del país. En esta zona 
norte, el acuífero se encuentra en contacto con 
formaciones poco permeables, por lo que fue 
considerado como un contorno impermeable. 

Las fronteras este y oeste del modelo pue- 
den ser asumidas como impermeables debido 
a que regionalmente existe un flujo del norte 
hacia el sur en estas zonas que no permite 
entradas ni salidas de agua subterránea al 
acuífero por dichas fronteras. 

El límite sur es considerado como una zona 
de descarga subterránea al mar, con cargas 
piezométricas conocidas en todos los nodos de 
este contorno. Se tomó para estos nodos una 
carga constante en el tiempo de 0.2 m, para así 
considerar el incremento de las cargas debido a 
la acción del dique sur construido a lo largo de 
la costa (ver Llanusa et al., 2004). 


Propiedades hidrogeológicas de la zona de 
estudio 


Se consideran como propiedades hidro- 
geológicas del acuífero: la conductividad 
hidráulica, K y la porosidad efectiva o coefi- 
ciente de almacenamiento S. Las propiedades 


Cuadro 1. Propiedades hidrogeológicas de los elementos del modelo. 


Material K (m/día) S Elementos Material K (m/día) 5) Elementos 
il 4.9 0.1 52, 53, 66, 68 20 789.0 0.2 34, 35 
2 9.9 0.1 59 228) 443.8 0.2 ZN, 2282, 209) 
3 29.6 0.1 60, 32 24 493.1 0.2 18, 19, 20 
4 49.3 0.1 27,30 225) 295.9 0.2 87 
5 74.0 0.2 50 27 1 534.2 0.2 5, 6,7, 8 
6 98.6 0.2 23, 24 28 1.0 0.1 67 
8 197.3 0.2 58, 31 29 3.0 0.1 64, 65, 69 
10 295:9) 0.2 47, 51, 28, 29, 26 30 4.9 0.2 11, 12, 49 
13 443.8 0.2 17, 44 31 2.6 0.1 55 
14 493.1 0.2 41, 42 S2 2.6 0.1 56 
15 542.5 0.2 ? EMO 16 33 23.0 0.1 61 
16 591.8 0.2 2 eS a 34 32m 0.1 62 
17 641.1 0.2 1 35 1.6 0.1 63 
19 VII 0.2 3) 36 14.0 0.1 57 
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hidrogeológicas de los elementos que confor- 
man la zona de estudio son obtenidas de la 
etapa de calibración dadas en 36 subzonas, 
columna Material, en el cuadro 1 (Llanusa et 
al., 1993). 


Niveles piezométricos del acuífero: estado 
inicial 


El estado inicial en el acuífero se obtuvo a 
partir de los valores medidos en los pozos 
de observación y su cota topográfica. Estos 
datos se tomaron de una serie de 19 pozos 
ubicados en la zona acuífera de Guúira- 
Quivicán y pertenecientes a la red informativa 
de niveles. 

Para definir el estado inicial en la modela- 
ción se tomaron los datos pertenecientes a la 
fecha (noviembre del año 1993), publicados 
en Hernández (2000); estos datos definen una 
carga media a lo largo de la costa de 0.57 m, 
producto de la presencia del dique sur. 


Alimentación del acuífero debido a la lluvia: 
recarga 


Batista (1995) hace un análisis de la lluvia 
media hiperanual desde el año 1964 hasta 
1993 con la información de 12 pluviómetros. 
De este análisis determinó la lluvia del 75% de 
probabilidad (1 113 mm), que prácticamente 
coincidió con la media obtenida en el año 
1975 (1 114 mm). Posteriormente tomó la 
distribución de lluvia de este año y realizó un 
balance hídrico diario del suelo de cobertura 
de la zona de estudio, empleando para esto 
el modelo matemático EVABAL. Finalmente 
obtuvo la distribución espacial de la recarga 
al acuífero debido a la lluvia, así como su 
cuantificación en volumen. Esta distribución 
quedó organizada en 12 grupos de triángulos, 
en los cuales la recarga mensual era diferente. 
Estos resultados se emplean en el presente 
trabajo, y se muestran en la figura 3 y en el 
cuadro 2. 

En el cuadro 2 aparecen los valores de 
recarga al acuífero expresados en forma de 
intensidad (mm/mes). Obsérvese que cada 


Figura 3. Distribución por grupos de la recarga subterránea 


al acuífero. 


columna representa un mes (AT) del año, que 
es el intervalo de tiempo con que se trabajará 
en este artículo, y cada fila está asociada con 
uno de los grupos de infiltración expuestos en 
la figura 3. 


Tratamiento de la explotación al acuífero 


Los volúmenes de explotación por nodos que 
serán empleados en las nuevas simulaciones 
fueron tomados de Batista (1995) y se muestran 
en el cuadro 3. 

En este cuadro se colocan primeramente las 
principales extracciones realizadas al acuífero, 
las cuales corresponden al acueducto y a la 
agricultura. La suma de estas extracciones 
asciende a 264.10 hm*/año, lo que representa 
un 90% de la extracción global, 295.20 hm?/ 
año. El resto está muy diseminado en otros 64 
nodos de la malla; en la última fila del cuadro 
3 se colocaron los acumulados de estos valores 
para este grupo de nodos. 


Análisis de la demanda de explotación global 
solicitada al acuífero 


Una forma de examinar si el acuífero puede 
satisfacer la explotación global solicitada es 
realizando una simulación con el criterio de 


Cabrera-Estupiñán y Dilla-Salvador, Modelo de administración de acuíferos: MADA 


Cuadro 2. Recarga al acuífero en lámina de agua acumulada (mm /mes) y por grupo de infiltración. 


Valores de recarga al acuífero por meses (mm/mes) 


Grupo AT, AT, AT, AT, AT, AT, AT, AT, AT, AT AT, AT, 
1 0 0 0 0 12.6 98.9 73.4 120.5 156.6 107.1 62.5 6.4 
2 0 0 0 0 0 22.5 108.4 109.8 110.6 32.8 1.2 0.1 
3 0 0 0 0 0 2.9 73 173 13.8 41.4 3.2 0.1 
4 0 0 0 0 0 0 0 0 0 2.8 42.8 2.7 
5 0 0 0 0 2.9 4.0 76.9 133.6 156.2 46.5 1.8 0.1 
6 0 0 0 0 14.8 56.7 28.7 62.4 72.2 623 2.4 0.1 
7 0 0 0 0 62.1 15.4 0.7 0 0 0 0 0 
8 0 0 0 0 25 10.3 0.5 0 33.7 61.8 2.3 0.1 
9 0 0 0 0 4.4 8.1 81.5 121.0 174.3 39.8 1.5 0.1 
10 0 0 0 0 2.0 60.0 47.8 164.8 229.5 64.9 4.5 0.2 
11 0 0 0 0 0 116.7 84.5 79.0 54.0 51.0 1.9 0.1 
12 0 0 0 0 39.8 59.9 76.2 249.3 151.7 44.6 1.7 0.1 


Cuadro 3. Nodos de explotación del modelo con sus volúmenes extraídos en cada mes. 


Nodos de explotación del acueducto, caudal en hm*/mes 


Nodos AT, AT, AT, AT, AT, AT, AT, AT, AT, AT AT. AT, 
45 0.16 0.13 0.11 0.12 0.28 0.29 0.27 0.35 0.00 0.46 0.07 0.09 
47 0.50 0.51 0.48 3.05 0.57 0.49 0.49 0.62 0.43 0.43 0.51 0.51 
7 0.86 0.69 0.85 0.90 0.99 0.74 0.91 0.80 0.73 0.80 0.89 0.83 
73 0.80 0.80 0.65 0.73 0.73 0.80 0.80 0.80 0.73 0.73 0.80 0.80 
85 0.40 0.38 0.45 0.48 0.44 0.45 0.36 0.36 0.36 0.36 0.45 0.45 
86 0.36 0.44 0.44 0.44 0.44 0.36 0.36 0.36 0.36 0.36 0.44 0.44 
87 0.36 0.58 0.52 0.51 0.46 0.44 0.29 0.29 0.36 0.36 0.44 0.44 
88 0.91 0.93 0.88 0.86 0.84 0.83 0.73 0.73 0.80 0.73 0.93 0.90 
101 0.40 0.49 0.47 0.43 0.40 0.46 0.36 0.36 0.36 0.36 0.44 0.36 
102 0.80 0.80 0.80 0.80 0.80 0.73 0.80 0.80 0.80 0.65 0.80 0.80 
103 0.91 0.86 1.00 1.01 0.88 0.84 0.88 0.81 0.80 0.78 0.78 0.80 
104 0.40 0.37 0.47 0.47 0.46 0.42 0.50 0.53 0.44 0.47 0.50 0.49 
119 0.73 0.80 0.80 0.80 0.80 0.73 0.65 0.80 0.80 0.80 0.80 0.80 
123 0.37 0.34 0.34 0.36 0.10 0.20 0.24 0.00 0.30 0.28 0.18 0.23 > 
124 0.36 0.33 0.40 0.39 0.35 0.23 0.23 0.00 0.07 0.16 0.26 0.16 5 
Nodos de explotación de la agricultura, caudal en hm?/mes » 
sl 0.08 0.45 0.09 0.03 0.01 0.03 0.67 SO 0.20 0.03 0.19 0.08 E 
4 0.85 0.86 0.45 0.00 0.32 0.50 0.16 0.86 0.09 1.01 0.46 0.45 + 
5 1.10 0.80 0.67 1.62 0.28 0.71 0.96 0.78 IZ 0.91 0.83 0.78 E 
13 0.85 0.55 0.90 0.17 0.68 0.00 0.64 0.87 0.44 0.05 0.00 0.09 E 
14 209 5.08 2.41 2.01 0.97 1.66 1.36 0.55 2.01 2.14 2.04 es £ 
15 1.19 0.32 0.14 0.39 0.65 0.04 0.76 0.16 0.59 0.24 0.89 0.12 = 
16 1.09 0.82 077 0.39 0.01 0.73 0.88 1.51 0.39 0.01 0.76 0.45 ? 
24 0.48 1.37 0.06 0.12 0.64 1.64 0.95 0.23 0.57 0.57 0.24 0.55 s 
25 3d lod 075 0.48 0.99 0.08 2.46 0.59 2139 1.90 1.58 29 3 
E 
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Cuadro 3 (continuación). Nodos de explotación del modelo con sus volúmenes extraídos en cada mes. 


Nodos de explotación de la agricultura, caudal en hm'/mes 
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26 0.73 0.24 7.96 7.77 3.32 1.08 0.64 0.50 0.13 0.48 0.88 0.46 
33 0.20 0.55 0.18 0.47 0.70 1.78 0.08 0.04 0.51 0.33 0.05 0.12 
36 0.62 0.78 0.16 0.23 0.09 0.17 0.47 0.25 0.95 1.02 1.02 0.73 
3 0.66 0.10 0.15 LT 077 1.24 0.08 0.04 0.02 0.05 0.11 1.49 
131 0.70 0.61 0.60 0.96 0.40 0.05 0.05 0.13 0.11 0.29 0.36 0.30 
135 0.16 0.19 0.14 0.48 0.09 0.10 0.06 0.05 0.06 0.06 0.06 0.05 
151 0.21 0.29 0.26 1.81 0.26 0.22 0.21 0.21 0.21 0.22 0.20 0.19 
152 0.24 0.24 0.24 0.49 0.03 0.00 0.01 0.04 0.03 0.00 0.00 0.02 
153 1.11 1.31 116 0.34 0.25 0.05 0.04 0.00 0.05 0.29 0.25 0.33 
154 0.64 0.74 193 1.03 0.21 0.09 0.07 0.01 0.06 0.06 0.07 OS 
158 1.36 0.92 0.87 1.84 0.38 0.06 0.04 0.00 0.00 0.06 0.37 0.67 
156 0.81 0.91 0.86 0.10 0.06 0.02 0.02 0.02 0.02 0.05 0.02 0.27 
Otras extracciones, caudal en hm*/mes 

Grupo 
de ano SS 2 4.05 2.56 20 2.48 2.05 2.68 2.07 191 2.42 

nodos 


que al final del año el volumen de agua en el 
acuífero sea similar al del inicio de la simula- 
ción, concepto equivalente al rendimiento 
seguro de la escuela norteamericana (Her- 
nández, 1991). 

En la figura 4a se muestra un gráfico de 
volúmenes en el acuífero en cada mes del 
año simulado con el modelo AQUIMPE. Es 
interesante observar que bajo el régimen de 
explotación y de recarga antes definidas, 
el acuífero al final del año queda con un 
volumen ligeramente menor al del inicio de la 
simulación. 

Un análisis más consistente pudiera ser 
el de repetir este año patrón varios años 
seguidos, para observar la tendencia de regu- 
lación del acuífero. Por ejemplo, se simularían 
cinco años de baja recarga al acuífero (75% de 
probabilidad), con una explotación constante 
(igual al PUA del año 1995). El gráfico de 
volúmenes de esta variante se muestra en la 
figura 4b. 

En la figura 4b se aprecia una tendencia 
de estabilidad del acuífero, lo que demuestra 
que el Plan de Uso del Agua de 295.20 hm?/ 
año, manejado por Batista (1995), tiene un 
gran valor práctico, al garantizar el equilibrio 


hidráulico en el acuífero bajo la lluvia del 75% 
probabilidad. 

En la figura 5 se muestra un mapa de 
niveles del agua subterránea correspondiente 
al tiempo cinco de las simulaciones realizadas. 
En la figura 4a se observa que es el tiempo en el 
primer año en que el acuífero acumuló menor 
volumen y por lo tanto se deben esperar en 
toda la región menores niveles piezométricos. 
De este mapa es importante resaltar que existe 
un conjunto de nodos, ubicados hacia la 
izquierda y en el centro, en donde las cargas 
hidráulicas están por debajo del nivel cero o 
nivel medio del mar. 

Este aspecto es extremadamente importante 
en acuíferos de este tipo, que están en contacto 
con el mar, debido a la posible contaminación 
de pozos, producto del avance de la intrusión 
salina provocada por el bombeo excesivo en la 
zona. 

Es por esta razón que en el presente trabajo 
se aplica un modelo de administración, para 
proponer una política de explotación que 
garantice una demanda global aceptable y que 
conserve los recursos hidráulicos. Además, 
debe estar distribuida de una forma eficaz en 
cuanto a garantizar que, en ningún pozo y en 
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Gráfico de volúmenes de agua en todo el acuífero Gráfico de volúmenes de agua en todo el acuífero 
para cada tiempo hm3 para cada tiempo 
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Tiempo (meses) Tiempo (meses) 
a) Un año b) Cinco años 


Figura 4. Gráfico de volúmenes por meses en todo el acuífero. 
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Figura 5. Mapa digital de isolíneas de niveles del agua subterránea en el tiempo cinco. 


Ciencias del Agua, vol. U, núm. 4, octubre-diciembre de 2011 


Tecnología y 


Ciencias del Agua, vol. Il, núm. 4, octubre-diciembre de 2011 


«o180/0u3aJ 


Cabrera-Estupiñán y Dilla-Salvador, Modelo de administración de acuíferos: MADA 


ningún tiempo, existan niveles por debajo del 
nivel medio del mar. 


Aplicación del modelo MADA en la 
optimización de la explotación del acuífero 


La primera variante analizada emplea un 
intervalo mensual, utilizando el mismo estado 
inicial y la recarga dados anteriormente. 
Se tendrá en cuenta una extracción global 
conocida, 295.20 hm? /año. 

Serán variables de decisión en el problema 
los caudales de los 36 nodos que se muestran 
en el cuadro 3 en cada uno de los tiempos y 
las cargas piezométricas en todos los nodos y 
tiempos. 

En el cuadro 4 se definen las extracciones 
mínimas de un grupo de nodos que tienen 
asociado caudal de extracción para acueductos, 
así como la máxima capacidad de extracción en 
ellos. 

Por otra parte, el cuadro 5 muestra los 
valores de extracción mínima y máxima de 
grupos de nodos vinculados con la explotación 
con fines agrícolas. Por ejemplo, la suma de los 
caudales extraídos en los nodos 3, 4 y 5 en cada 
tiempo está entre 0.01 y 5 hm*/mes. 


Cuadro 4. Límites inferior y superior de caudal de los 


nodos asociados con los acueductos. 


Nodos Extracción mínima | Extracción máxima 
(hm*/mes) (hm? /mes) 

45 0.1 0.5 

47 0.4 510 

71 0.7 1 

73 0.65 0.8 

85 0.36 0.50 

86 0.36 0.5 

87 0.36 0.6 

88 0.72 1 

101 0.36 0.5 

102 0.8 0.8 

103 0.8 1.1 

104 0.36 0.6 

119 0.8 0.8 

125 0.1 0.4 

124 0.1 0.4 


Cuadro 5. Límites inferior y superior de caudal para grupos 
de nodos asociados con la agricultura. 


Extracción mín. | Extracción máx. 
Grupo de dl 1 
lO noda grupa 
hm'/mes hm'/mes 

3,45 0.01 5 
13, 14, 15, 16 0.05 6 
24, 25, 26 0.05 5 
SOMOS 0.1 3 
131, 135 0.05 2 
151 al 156 0.02 8 


Cuadro 6. Límites inferior y superior de los niveles 
en metros sobre el nivel medio del mar (msnmm) en un 
grupo de nodos. 


Nodos Nivel mínimo Nivel máximo 
(msnmm) (msnmm) 

5 0.50 00 

28 0.50 00 

44 0.40 00 

46 0.40 00 

48 1.30 00 

67 0.35 00 

99 0.20 00 
105 1.00 00 
107 4.00 00 
lo) 0.10 00 
132 1.00 00 


Los límites de los niveles para todos los 
nodos en todos los tiempos, exceptuando los 
que se encuentran en el cuadro 6 coincidentes 
con pozos de observación, están definidos entre 
0 e infinito (00). 

Para cada uno de los intervalos de tiempos 
(meses) se define un valor de demanda acorde 
con los datos de extracciones manejados 
del PUA de 295.20 hm*/año, estos valores 
organizados por meses, en hm*/mes, son lo 
siguientes: 28.0, 28.0, 28.0, 34.0, 20.0, 18.0, 18.5, 
18.0, 16.5, 17.5, 18.6, 19.0. 

Obsérvese que la suma de ellos da un valor 
de 264.10 hm? /año, el cual es menor que el del 
PUA y esto tiene sentido, ya que la diferencia 
es de alrededor del 10% y corresponde con un 
caudal que no será optimizado, o sea que es 
conocido e impuesto, última fila del cuadro 3. 
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La función objetivo utilizada en este 
estudio es la que responde a la ecuación (11), o 
sea, implica maximizar la suma de los caudales 
en los 36 nodos de extracción, en todos los 
tiempos. 

En esta primera optimización, las res- 
tricciones de demanda se definieron como 
menor o igual que. Esto garantiza, por ejemplo, 
que en el primer tiempo, la suma de todos 
los caudales en los 36 nodos de optimización 
sea como máximo igual a 28.0 hm?/mes, y 
como la función objetivo implica maximizar el 
caudal, es de esperar que ese sea el resultado. 
En este caso, es interesante el resultado de la 
distribución que propone el modelo, que no 
tiene que coincidir con la forma de explotación 
real, dada a priori al acuífero, expuesta en el 
cuadro 3. 

Los resultados de los caudales, una vez 
ejecutado el modelo MADA, se muestran en 
el cuadro 7 y los volúmenes obtenidos en el 
acuífero para cada tiempo en la figura 6a. 

Es importante recordar que con esta 
distribución de caudales se cumplen todas las 


restricciones definidas anteriormente, lográn- 
dose que no existan niveles negativos en 
ninguna zona del acuífero. 

Se realizó otra optimización con las mismas 
características que la anterior, pero en este 
caso sin tomar en cuenta las restricciones 
de demanda, es decir, que el modelo podría 
decidir extraer toda la máxima cantidad de 
agua posible cada mes, garantizando los niveles 
piezométricos en los límites antes definidos. 
Los resultados gráficos en volúmenes de esta 
simulación se muestran en la figura 6b. 

Como puede observarse en la figura 6b, los 
volúmenes acumulados son menores y esto 
se debe a una fuerte extracción, que ascendió 
a unos 339.28 hm?*/año, lo que representa un 
128.5% del volumen necesario (264.10 hm?/ 
año) para garantizar el PUA. Esta excesiva 
explotación lleva al acuífero a terminar en 
una situación crítica en los niveles al fin del 
año. 

Este resultado obtenido tiene un gran 
valor práctico, ya que brinda una idea del 
máximo caudal anual posible a extraer bajo las 


Gráfico de volúmenes de agua en todo el acuífero 
para cada tiempo 


Gráfico de volúmenes de agua en todo el acuífero 
para cada tiempo 


Tiempo (meses) 


Tiempo (meses) 


b) 


Figura 6. a) Volúmenes en el acuífero obtenidos en la optimización con restricciones en la demanda; 


b) volúmenes en el acuífero obtenidos en la optimización sin restricciones en la demanda. 
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Cuadro 7. Caudales optimizados para cada tiempo en los 36 nodos del modelo de administración. 


Nodos de explotación del acueducto, caudal en hm'/mes 


Nodos AT, AT, AT, AT, AT, AT, AT, AT, AT, AT, AT, AT, 
45 0.07 0.07 0.04 0.42 0.26 0.26 0.25 0.22 0.21 0.23 0.24 0.26 
47 0.08 0.48 0.01 2.71 0.33 1.32 1.17 0.88 0.65 0.74 0.84 0.88 
Al 0.23 0.13 0.32 0.29 0.18 0.27 0.31 0.28 0.27 0.31 0.35 0.39 
73 0.56 0.73 0.73 0.76 0.60 0.50 0.45 0.35 0.31 0.33 0.36 0.37 
85 0.21 0.18 0.24 0.34 0.18 0.23 0.23 0.21 0.20 0.22 0.23 0.25 
86 0.22 0.37 0.38 0.44 0.27 0.28 0.27 0.23 0.22 0.23 0.24 0.25 
87 0.43 0.54 0.53 (057 0.44 0.39 0.36 0.30 0.27 0.28 0.30 0.31 
88 0.96 0.96 0.96 0.97 0.86 0.76 0.70 0.54 0.46 0.47 0.48 0.49 
101 0.19 0.15 0.21 0.16 0.15 0.23 0.24 0.21 0.20 0.22 0.23 0.24 
102 0.16 0.67 0.72 0.71 0.45 0.46 0.43 0.33 0.30 0.33 0.35 0.37 
103 0.87 1.03 1.03 1.06 0.84 0.73 0.64 0.45 0.38 0.41 0.44 0.46 
104 0.55 0.55 0.55 0.57 0.48 0.43 0.40 0.33 0.30 0.31 0.32 0.32 
119 0.65 0.63 0.68 0.70 0.46 0.46 0.43 0.34 0.31 0.33 0.36 0.37 
123 0.38 0.38 0.38 0.38 0.31 0.30 0.28 0.24 0.22 0.22 0.22 0.21 
124 0.38 0.37 0.37 0.38 0.33 0.30 0.28 0.24 0.22 0.22 0.22 0.21 

Nodos de explotación de la agricultura, caudal en hm*/mes 

5) 0.67 0.35 0.22 0,12 0.07 0.09 0.09 0.08 0,11 0.21 0.19 0.07 
0.67 0.23 0.26 0.94 0.11 0.14 0.31 0.39 0.42 0.48 0.57 0.68 

5 0.05 0.16 0.46 0.18 0.02 0.03 0.06 0.34 0.37 0.33 0.23 0.14 
13 2.50 1.52 0.36 0.93 0.10 0.13 0.08 0:25 0.30 0.38 0.44 0.59 
14 0.16 0.24 0.13 0.91 0.05 0.07 0.24 0.44 0.42 0.45 0.41 0.51 
15 0.27 0.46 0.24 2.36 0.09 0.15 0.55 0.76 0.60 0.58 0.58 0.80 
16 1079 114122 2.08 1.14 0.45 0.42 0.65 0.66 0.58 0.64 0.83 0.67 
24 0.95 0.20 1.25 0.51 0.12 0.27 0.31 0.43 0.45 0.55 0.70 0.61 
5 0.44 0.19 MIT 2.90 0.49 0.57 1152 1.05 0.83 0.88 1.10 1.20 
26 1.81 3.63 1.41 0.80 0.53 0.60 0.76 0.98 0.89 0.96 0.95 0.74 
35 0.04 0.07 0.04 0.04 0.12 0.36 0.40 0.50 0.50 0.57 0.61 0.63 
36 0.08 0.13 0.06 0.18 0.60 0.82 0.94 0.74 0.65 0.68 0.71 0.73 
SY 2.80 2.63 279 2 1.90 112 0.99 0.80 0.69 0.69 0.69 0.66 
131 0.07 0.07 0.04 0.42 0.26 0.26 0.25 0.22 0.21 0.23 0.24 0.26 
185 0.08 0.48 0.01 JA 0.33 11.592 Li 0.88 0.65 0.74 0.84 0.88 
151 0.23 0.13 0.32 0.29 0.18 0.27 0.31 0.28 0.27 0.31 0.35 0.39 
152 0.56 0.73 0.73 0.76 0.60 0.50 0.45 0.35 0.31 0.33 0.36 0.37 
153 0.21 0.18 0.24 0.34 0.18 0.23 0.23 0.21 0.20 0.22 0.23 0.25 
154 0.22 0.37 0.38 0.44 0.27 0.28 0.27 025 122 0.23 0.24 0.25 
155 0.43 0.54 0.53 0.57 0.44 0.39 0.36 0.30 0.27 0.28 0.30 0.31 
156 0.96 0.96 0.96 0.97 0.86 0.76 0.70 0.54 0.46 0.47 0.48 0.49 


condiciones de un año seco (75% de probabi- 
lidad de lluvia) y garantizando niveles piezo- 
métricos por encima del nivel medio del mar. 
Finalmente, se realizó una optimización 
de la explotación, repitiendo el mismo 
año hidrológico cinco veces, por lo que se 


tomaron en cuenta en total sesenta intervalos 
de tiempo. Lo que se busca ahora con este 
análisis es apreciar el comportamiento 
del acuífero bajo la acción de la demanda 
solicitada en un periodo seco prolongado, en 


este caso de cinco años. 
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Los resultados obtenidos de la simulación 
en cuanto al cumplimiento de la demanda 
exigida cada mes en cada uno de los cinco años 
se presentan en el cuadro 8. En dicho cuadro 
se observa que en el primer año se garantiza 
la demanda planteada, sin embargo en los 
meses tres y cuatro del segundo año, no se 
puede entregar la demanda y es lógico que 
sea precisamente en estos meses que comience 
la crisis, ya que en los primeros cuatro meses 
del año no hay recarga al acuífero (ver cuadro 
2), sin embargo, se le pide una extracción 
considerable, sobre todo en el mes cuatro, que 
asciende a 34.0 hm" y este efecto puede ser 
soportado el primer año, pero posteriormente 
todos los años presentan déficit. Los meses con 
déficit en cada año pueden observarse en el 
cuadro 8, en las celdas marcadas. 

Con estos resultados se demuestra que el 
acuífero no puede cumplir con las demandas 
mensuales planteadas, por lo tanto éstas 
deberán disminuir, de forma tal que se pueda 
encontrar una distribución máxima de las 
demandas que garantice la no ocurrencia de 
años con déficit. En este caso, la distribución 
encontrada mediante un ensayo de prueba y 
error es la siguiente 28.0, 23.0, 20.0, 20.0, 19.0, 
18.0, 18.5, 18.0, 16.5, 17.5, 18.6, 19.0. 


La suma de estos valores es de 236.10 hm?/ 
año. La diferencia entre el PUA y este valor es 
de 28.0 hm*/año. Pudiera pensarse en variar 
la distribución en función de los objetivos 
que se quieran garantizar, sin embargo, con 
esta distribución de demandas por meses que 
se Obtuvo, se garantiza la entrega para las 
condiciones exigidas. 


Análisis de dimensiones y tiempos de 
ejecución con el modelo MADA 


Un aspecto interesante en la última ejecución 
del modelo MADA está referido a la dimensión 
de la matriz A,, que contiene al modelo 
AQUIMPE en el algoritmo de programación 
lineal. Para este ejemplo, dicha matriz necesitó 
de 9 480 filas y 11 640 columnas para un total 
de 110 347 200 elementos. El programa está 
preparado a tal efecto para trabajar con matrices 
de tipo “sparse”, con las cuales se reduce el 
espacio de almacenamiento de los datos, ya 
que se almacenan sólo los valores diferentes de 
cero de las matrices. 

A pesar de las dimensiones del problema 
solucionado, el método solamente empleó 
unos dos minutos, con un total de 49 
iteraciones en un ordenador personal con 


Cuadro 8. Valores de los volúmenes de agua extraídos cada mes en cada uno de los cinco años simulados. 


Extracciones acumuladas (hm*/año) 
Mess Demandas prefijadas 
(hm*/año) Año 1 Año 2 Año 3 Año 4 Año5 
1 28 28 28 28 26.23 26.85 
2 28 28 28 19.36 18.49 18.13 
3 28 28 2581 20.39 20.37 95) 
4 34 8979 20.7 sion 13.89 529] 
5 20 20 20 IIS 18.99 20 
6 18 18 18 18 18 18 
7 18.5 18.5 18.5 18.5 18.5 18.5 
8 18 18 18 18 18 18 
9 16.5 16.5 16.5 16.5 16.5 16.5 
10 AS MAS 11745 1745 11745 17/05 
1 18.6 18.6 18.6 18.6 18.6 18.6 
1 19 19 19 19 19 19 
Total 264.1 263.89 248.11 226.91 224.07 224.32 
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microprocesador Intel Core 2 Duo a 2.40 GHZ 
y 2 GB de memoria RAM, con lo cual se puede 
concluir que cuando el problema a resolver 
tiene grandes dimensiones —si bien el método 
de programación lineal necesita almacenar 
mucha información— el resultado se obtiene en 
un tiempo muy corto y con la garantía de que 
si el problema está bien planteado la solución 


será la correcta. 


Conclusiones 


1. Se realiza la concepción teórica general 
de un modelo de optimización para 
la administración de la explotación de 
acuíferos  hidráulicamente confinados, 
llamado MADA. La originalidad del 
modelo propuesto se basa en introducir 
las ecuaciones, que emplea el modelo de 
simulación de acuíferos AQUIMPE, como 
restricciones de igualdad, al emplear la 
programación lineal como técnica de 
optimización. Además, se tienen en cuenta 
un conjunto de restricciones a las variables 
de decisión que enriquecen el valor práctico 
de este modelo. 

2. Queda demostrada la factibilidad del 
modelo propuesto MADA, al llevarse a cabo 
una evaluación crítica de la explotación 
en el modelo regional de la Cuenca Sur 
de La Habana. En este caso de estudio se 
demostró que los resultados que brinda el 
modelo MADA para la distribución de la 
explotación que garantiza el Plan de Uso del 
Agua, evalúa las posibilidades de intrusión 
salina en el acuífero, o sea, se obtiene una 
distribución de caudales que respeta los 
valores piezométricos regionales críticos en 
el acuífero. 

3. Se comprueba que el nivel de extracciones 
correspondiente al Plan de Uso del Agua 
puede verse afectado por el efecto de varios 
años secos continuados. En este sentido se 


realiza una optimización de las extracciones 
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para cinco años secos consecutivos. Los 
resultados muestran años con meses en los 
que no se cumple la demanda (meses con 
déficit). Además se muestra la flexibilidad 
del modelo para encontrar soluciones 
a estos problemas, dependiendo de los 
objetivos que persiga el modelador. 
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Abstract 


CABRERA-ESTUPINÁN, E. és DILLA-SALVADOR, F. Groundwater management model: 
MADA. Water Technology and Sciences (in Spanish). Vol. II, No. 4, October-December, 
2011, pp. 5-24. 


In this paper a mathematical model (MADA) is proposed for the optimal management of aquifer 
exploitation. This model is based on the inclusion of a simulation model for groundwater flow 
(AQUIMPE) using a linear programming algorithm. The theoretical concept behind MADA 
is described, making a special distinction for the inclusion of AQUIMPE in the proposed 
method. Finally, MADA is used to evaluate the exploitation of the South Havana Basin 
according to a Water Use Plan (WUP), where it is demonstrated that the distribution of the 
exploitation used in the WUP does not take into account possible saline intrusion into the 
aquifer; a solution to this effect is provided. Also evaluated is the ability of the aquifer to satisfy 
demand based on a yearly pattern of75% probability of rain during a period of five consecutive 
years with low rainfall. 


Keywords: optimization, ground-water modeling, AQUIMPE, MADA, linear programming, 
South Havana Basin. 
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